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ABSTRACT 

A  new  modeling  framework  is  introduced  for  the  analytical 
study  of  medium  access  control  (MAC)  protocols  operating 
in  multihop  ad  hoc  networks.  The  model  takes  into  ac¬ 
count  the  effect  of  physical-layer  parameters  on  the  success 
of  transmissions,  the  MAC  protocol  on  the  likelihood  that 
nodes  can  access  the  channnel,  and  the  connectivity  of  nodes 
in  the  network.  A  key  feature  of  the  model  is  that  nodes  can 
be  modeled  individually,  i.e.,  it  allows  a  per-node  setup  of 
many  layer-specific  parameters.  Moreover,  no  spatial  prob¬ 
ability  distribution  or  a  particular  arrangement  of  nodes  is 
assumed;  the  model  allows  the  computation  of  individual 
(per-node)  performance  metrics  for  any  given  network  topol¬ 
ogy  and  radio  channel  model.  To  show  the  applicability  of 
the  modeling  framework,  we  model  multihop  ad  hoc  net¬ 
works  using  the  IEEE  802.11  distributed  coordination  func¬ 
tion  and  validate  the  results  from  the  model  with  discrete- 
event  simulations  in  Qualnet.  The  results  show  that  our 
model  predicts  results  that  are  very  close  to  those  attained 
by  simulations,  and  requires  seconds  to  complete  compared 
to  several  hours  of  simulation  time. 
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1.  INTRODUCTION 

The  medium  access  control  (MAC)  protocol  of  a  computer 
network  enables  nodes  to  determine  their  right  to  access  the 
available  channel(s),  while  attempting  to  enforce  fair  and 
efficient  usage  of  the  channel(s).  Establishing  such  access 
rights  is  far  more  difficult  in  an  ad  hoc  network  than  in  a 
wired  long  haul  network  or  a  local  area  network  (LAN),  be¬ 
cause  the  radio  channels  of  an  ad  hoc  network  are  broadcast 
in  nature  and  radio  connectivity  is  such  that  the  topology  of 
an  ad  hoc  network  is  not  as  clearly  defined  as  with  point-to- 
point  wire  networks.  Unlike  wired  channels  that  are  fairly 
stationary  and  predictable,  radio  channels  are  extremely 
random,  and  connectivity  between  two  nodes  depends  on 
many  factors,  such  as  the  radio  frequency  in  use,  power  of 
the  transmitters,  terrain,  antenna  type,  transmitter/receiver 
distance,  multipath  fading,  and  the  like.  Furthermore,  the 
quality  of  a  radio  link  depends  on  the  transmission  activ¬ 
ity  of  all  other  nodes  in  the  entire  system,  whose  aggregate 
signal  powers  can  severely  degrade  the  signal-to-noise  ratio 
(SNR)  at  a  particular  receiver  and,  consequently,  compro¬ 
mise  the  successful  reception  of  any  on-going  packet  trans¬ 
mission.  Last,  but  not  least,  thermal  and  other  background 
noise  sources  can  also  contribute  to  the  failure  of  a  packet 
reception  at  a  given  receiver. 

Hence,  analyzing  the  performance  of  a  MAC  protocol  in 
an  ad  hoc  network  must  consider  the  interactions  between 
the  physical  (PHY)  and  the  MAC  layers.  In  fact,  a  cross¬ 
layer  perspective  to  both  analysis  and  design  of  multiac¬ 
cess  communications  has  been  brought  to  attention  with  re¬ 
cent  advances  in  wireless  communications.  As  pointed  out 
by  Gallager  [1]  and  later  by  Ephremides  and  Hajeck  [2], 
the  fundamental  challenge  lies  in  the  choice  of  a  proper 
model  that  interfaces  the  physical  layer  and  network  lay¬ 
ers.  Unfortunately,  the  bulk  of  the  published  work  on  ad 
hoc  networks  has  focused  on  the  modeling  of  medium  ac¬ 
cess  schemes  where  ideal  channel  conditions  and  unrealistic 
assumptions  are  made.  Under  the  argument  of  separating 
the  issues  strictly  related  to  the  protocol  operation  from  the 
issues  intrinsically  related  to  the  physical  layer,  few  studies 
have  attempted  to  incorporate  physical  layer  aspects  directly 
into  the  behavior  of  the  protocol,  i.e.,  explicitly  modeling  the 
impact  of  the  physical  layer  on  the  dynamics  of  the  MAC 
protocol.  Perhaps  more  importantly,  no  attempt  has  ever 
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been  made  to  include  not  only  the  physical-layer  aspects 
explicitly  into  the  model,  but  also  to  include  the  interde¬ 
pendencies  among  nodes  under  a  radio-based  topology  in  a 
multihop  ad  hoc  network.  Although  convenient  for  analyti¬ 
cal  modeling,  it  is  not  true  that  only  the  transmissions  from 
“one  hop”  neighbors  of  a  node  can  cause  interference  at  the 
node. 

Because  of  the  limitations  of  existing  analytical  models 
for  ad  hoc  networks,  many  researchers  have  opted  to  study 
the  impact  of  physical  layer  on  the  dynamics  of  MAC  proto¬ 
cols  via  discrete-event  simulations.  However,  even  with  such 
simulation  packages  as  Qualnet,  which  are  designed  to  scale 
with  the  number  of  nodes  in  the  MANET,  obtaining  sim¬ 
ulation  results  that  are  statistically  meaningful  (by  using 
many  seeds)  requires  hours  of  simulation  time  for  scenar¬ 
ios  corresponding  to  just  a  few  minutes  of  simulated  time 
in  MANETs  with  hundreds  of  nodes  with  a  given  choice 
of  physical-layer  parameters.  This  is  clearly  not  a  promis¬ 
ing  approach  for  researchers  to  gain  insight  on  the  impact  of 
multiple  physical-layer  parameters  on  the  operation  of  MAC 
protocols! 

This  paper  introduces  a  new  modeling  framework  for  any 
MAC  protocol  operating  in  multihop  ad  hoc  networks  that 
focuses  on  the  interoperability  between  the  PHY  and  the 
MAC  layers.  In  this  sense,  the  model  builds  up  on  their  ulti¬ 
mate  functionalities,  which  are  those  of  providing  a  schedul¬ 
ing  discipline  for  nodes  to  have  access  to  the  channel,  and 
guaranteeing  that  a  transmitted  frame  is  correctly  received. 
To  account  for  the  effects  of  both  cross-layer  interactions 
and  the  interference  among  all  nodes,  a  novel  linear  model 
is  introduced  with  which  topology  and  physical  layer  aspects 
are  naturally  incorporated  in  what  we  define  as  interference 
matrices.  A  key  feature  of  the  model,  which  is  introduced 
in  Section  3,  is  that  nodes  can  be  modeled  individually, 

1. e.,  it  allows  a  per-node  setup  of  many  layer-specific  pa¬ 
rameters.  Moreover,  no  spatial  probability  distribution  or 
pre- arrangement  of  nodes  is  assumed.  On  the  contrary,  the 
model  allows  the  computation  of  individual  (per-node)  per¬ 
formance  metrics  for  any  given  network  topology  and  radio 
channel  model. 

Section  2  presents  an  overview  of  past  work  in  the  model¬ 
ing  and  performance  analysis  of  ad  hoc  networks.  Section  4 
illustrates  the  applicability  of  the  new  modeling  framework 
by  applying  it  to  the  well-known  IEEE  802.11  distributed 
coordination  function  (DCF).  In  doing  so,  new  analytical 
results  are  derived  that  extend  our  prior  work  on  fully- 
connected  networks  [3],  and  the  work  by  Bianchi  [4].  Sec¬ 
tion  5  validates  the  results  obtained  with  the  new  modeling 
framework  using  the  Qualnet  v3.5  simulation  package  [5], 
and  demonstrates  that  the  new  framework  produces  results 
that  correlate  quite  nicely  with  the  results  obtained  using 
simulations,  and  accomplishes  this  in  seconds  or  fractions 
of  a  second,  representing  five  or  six  orders  of  magnitude  in 
time  savings  with  respect  to  the  time  needed  to  run  the 
simulations.  Following  that,  Section  6  presents  an  exam¬ 
ple  of  how  the  new  modeling  framework  can  be  used  for  a 
fast  and  efficient  performance  evaluation  of  multihop  ad  hoc 
networks.  Finally,  Section  7  presents  our  final  remarks  and 
conclusions. 

2.  RELATED  WORK 

The  bulk  of  the  analytical  modeling  of  wireless  ad  hoc 
networks  has  concentrated  on  the  analysis  of  MAC  pro¬ 


tocols  in  fully-connected  segments  of  networks  (e.g.,  satel¬ 
lite  networks,  cellular  networks,  or  single-hop  wireless  LANs 
(WLANs)),  because  they  are  simpler  to  analyze  than  mul- 
tihop  networks.  The  majority  of  this  work  has  followed  the 
formalism  and  assumptions  introduced  by  Abramson  [6,  7] 
for  the  analysis  of  the  ALOHA  protocol,  and  by  Tobagi  and 
Kleinrock  [8,  9]  for  the  analysis  of  the  carrier  sense  multiple 
access  (CSMA)  protocol.  The  model  typically  adopted  as¬ 
sumes  that  all  nodes  have  infinite  buffers  and  transmissions 
are  scheduled  according  to  independent  Poisson  point  pro¬ 
cesses.  This  implies  that  packets  which  were  either  inhibited 
from  being  transmitted  or  were  unsuccessfully  transmitted 
are  rescheduled  after  a  “sufficiently  long”  randomized  time 
out  to  preserve  the  Poisson  property  (i.e.,  no  correlation  be¬ 
tween  new  packet  arrivals  and  their  rescheduling).  Packet 
lengths  are  exponentially  distributed  and  are  independently 
generated  at  each  transmission  attempt  (including  retrans¬ 
missions).  In  many  cases,  acknowledgments  are  assumed 
to  happen  instantaneously  or,  in  cases  where  propagation 
delay  is  taken  into  account,  acknowledgment  traffic  is  sim¬ 
ply  ignored,  and  periods  of  collisions  are  restricted  to  the 
propagation  time,  after  which  all  other  nodes  are  able  to 
perceive  any  activity  in  the  channel  (through  the  single-hop 
and  perfect-channel  assumptions).  Regarding  the  quality  of 
the  radio  links,  they  are  generally  considered  error  free,  and 
the  event  of  unsuccessful  transmission  is  restricted  to  packet 
collisions  at  the  receiver.  Examples  where  such  assumptions 
have  been  made  include  [10],  [11],  [12],  [13],  [14],  [15],  [16], 
and  [17]. 

Other  works  consider  physical-layer  aspects  more  explic¬ 
itly  within  the  context  of  single-hop  scenarios.  Raychauduri 
[18]  analyzed  slotted  ALOHA  with  code  division;  Grone- 
meyer  and  Davis  [19]  considered  spread-spectrum  slotted 
ALOHA  with  capture  due  to  time  of  arrival.  Musser  and 
Daigle  [20]  derived  the  throughput  of  pure  ALOHA  with 
code  division.  Pursley  [21]  studied  the  throughput  of  fre¬ 
quency-hopped  spread-spectrum  communications  for  packet 
radio  networks.  In  other  cases,  the  error-free  link  assump¬ 
tion  was  relaxed  and  multipath  fading  channels  where  con¬ 
sidered  while  preserving  other  original  assumptions  (e.g., 
Poisson  scheduling).  This  is  the  case  in  the  works  by  Arn- 
bak  and  Blitterswijk,  who  studied  the  capacity  of  slotted 
ALOHA  in  Rayleigh-fading  channels  [22]. 

More  recently,  with  the  advent  of  the  IEEE  802.11  stan¬ 
dard  for  WLANs,  a  significant  amount  of  work  has  been 
carried  out  to  provide  an  analytical  model  for  its  operation. 
Unfortunately,  the  vast  majority  of  this  effort  has  considered 
only  single-hop  networks  under  ideal  channel  conditions  [3] , 
[4],  [23],  [24],  [25].  As  far  as  imperfect  channel  conditions 
is  concerned,  Hadzi-Velkov  and  Spasenovski  [26,  27]  have 
investigated  the  impact  of  capture  on  the  capacity  of  the 
IEEE  802.11  in  both  Rayleigh-  and  Rician-fading  channels. 
However,  no  provision  was  made  to  consider  a  multihop  ad 
hoc  network  and  the  interdependencies  among  the  nodes. 

Gitman  [28]  published  what  is  arguably  the  first  paper 
that  actually  dealt  with  a  multihop  system.  Gitman  con¬ 
sidered  a  two-hop  centralized  network  consisting  of  a  large 
number  of  terminals  communicating  with  a  single  station  via 
some  repeaters  located  around  the  station.  Subsequently, 
Tobagi  [29,  30]  considered  the  same  topology  to  compute 
the  network  capacity  and  throughput-delay  characteristics 
of  both  slotted  ALOHA  and  CSMA.  However,  it  was  not 
until  the  work  by  Boorstyn  et  al.  [31]  that  a  methodology  for 


the  steady-state  throughput  analysis  of  a  multihop  packet 
radio  network  was  introduced.  Based  on  assumptions  like 
independent  Poisson  scheduling  at  each  node,  exponentially- 
distributed  packet  lengths,  zero  propagation  delay,  and  in¬ 
stantaneous  acknowledgments,  Boorstyn  et  al.  were  able  to 
represent  a  CSMA  multihop  network  as  a  continuous-time 
Markov  chain,  with  the  state  at  each  time  being  the  set  of 
transmitting  nodes.  This  analysis  lead  to  a  product-form 
solution  and  an  iterative  procedure  was  used  to  obtain  the 
scheduling  rates  corresponding  to  given  desired  link  traffic 
rates.  The  complexity  of  the  algorithm  by  Boorstyn  et  al., 
although  exponential  in  general,  grows  quadratically  or  cu- 
bically  with  the  number  of  links  for  most  networks  on  the 
order  of  100  nodes.  Tobagi  and  Brazio  [32]  considered  the 
same  model  and  observed  that  it  is  applicable  to  ALOHA 
and  C-BTMA  (a  variant  of  BTMA).  They  also  observed  that 
not  all  schemes  could  be  modeled  by  simply  tracking  the  set 
of  transmitting  nodes.  Shepard  [33]  considered  fundamen¬ 
tal  physical-layer  aspects  in  the  modeling  of  large,  dense 
packet-radio  networks.  Although  Shepard’s  model  consid¬ 
ered  aspects  closer  to  the  underlying  physics  of  radio  com¬ 
munications,  it  did  not  target  the  modeling  of  the  interaction 
among  the  nodes  and  their  interdependencies. 

Chhaya  and  Gupta  provided  one  of  the  first  analytical 
models  of  the  IEEE  802.11  DCF  that  considered  a  multi¬ 
hop  scenario  and  where  both  capture  and  hidden  terminals 
were  taken  into  account  [34].  Like  prior  work,  however,  all 
nodes  were  assumed  to  transmit  independently  according 
to  some  Poisson  point  process.  Moreover,  in  their  numer¬ 
ical  results,  they  considered  that  all  individual  scheduling 
processes  offered  the  same  load,  regardless  of  network  topol¬ 
ogy,  therefore  ignoring  the  many  interdependencies  among 
nodes  and  their  impact  on  individual  scheduling  rates.  Sub¬ 
sequently,  Wang  and  Garcia-Luna-Aceves  [35]  provided  a 
model  for  the  saturation  throughput  of  collision  avoidance 
protocols  that  included  the  IEEE  802.11  DCF.  Their  model¬ 
ing  approach  is  based  on  the  work  by  Wu  and  Varshney  [36] , 
in  which  nodes  are  spatially  distributed  according  to  a  two- 
dimensional  Poisson  distribution,  and  two  Markov  models 
are  used  to  represent  the  channel  around  a  node  and  its 
activity.  To  obtain  the  Markov  models,  the  channel  is  mod¬ 
eled  as  a  circular  region  in  which  nodes  within  the  region 
can  communicate  with  each  other,  while  weak  interactions 
are  assumed  with  nodes  outside  this  region. 

Another  important  venue  of  research  emerged  in  the  last 
few  years  that  focuses  on  the  problem  of  network  capacity. 
The  main  objective  of  this  line  of  work  consists  of  finding 
fundamental  limits  on  achievable  communication  rates  in 
wireless  networks.  In  such  a  problem  formulation,  a  set  of 
rates  between  source-destination  pairs  is  called  achievable  if 
there  exists  a  network  control  policy  that  guarantee  those 
rates.  The  closure  of  the  set  of  achievable  rates  is  the  ca¬ 
pacity  region  of  the  network.  In  their  seminal  work,  Gupta 
and  Kumar  [37]  considered  a  joint  optimization  of  trans¬ 
mission  powers  and  schedules,  showing  that  the  maximum 
per-node  throughput  scales  proportionally  to  1/ ^ Jn ,  where 
n  is  the  number  of  nodes.  Subsequent  works  have  studied 
network  capacity  from  various  viewpoints  (e.g.,  [38,  39,  40, 
41,  42,  43])  while  others  considered  the  joint  optimization  of 
resource  allocation  and  scheduling  [44,  45],  or  the  optimal 
cross-layer  design  of  PHY  and  MAC  layers  [46,  47] .  Despite 
the  undeniable  importance  of  such  studies,  a  gap  still  re¬ 
mains  on  the  modeling  of  wireless  ad  hoc  networks  under 


specific  (optimally  designed  or  not)  MAC  and  PHY  layers, 
in  a  way  that  the  impact  of  their  interactions  and  the  in¬ 
terdependence  among  the  nodes  in  a  multihop  environment 
are  all  taken  into  account  in  the  performance  evaluation  of 
each  node. 

3.  THE  ANALYTICAL  MODEL 

Our  analytical  model  focuses  on  the  essential  aspects  of 
the  PHY  and  MAC  layers,  as  well  as  their  interaction.  On 
the  one  hand,  the  MAC  layer  provides  a  scheduling  discipline 
for  nodes  to  access  the  shared  channel(s)  of  the  network,  and 
this  discipline  renders  probabilities  that  nodes  will  attempt 
to  transmit.  On  the  other  hand,  the  PHY  layer  encodes  in¬ 
formation  attempting  to  ensure  that  the  transmitted  frames 
are  received  correctly;  the  likelihood  with  which  a  transmis¬ 
sion  is  successful  depends  on  how  well  the  signaling  used  de¬ 
fends  against  channel  impairments  and  interference  from  any 
source.  Clearly,  the  dynamics  of  the  MAC  layer  is  tightly 
connected  to  the  dynamics  of  the  PHY  layer!  This  inter¬ 
action  between  PHY  and  MAC  layers  is  a  function  of  the 
connectivity  of  nodes  in  the  ad  hoc  network,  which  is  not  a 
deterministic  boolean  function,  but  rather  a  function  of  the 
location  of  nodes,  their  transmission  and  sensing  ranges,  and 
other  factors  that  together  determine  the  topology  of  the  net¬ 
work.  Accordingly,  the  performance  experienced  by  a  node 
is  a  complex  function  of  the  signals  used  at  the  PHY  layer, 
the  scheduling  established  at  the  MAC  layer,  and  the  topol¬ 
ogy  of  the  network. 

In  the  sequel,  we  build  our  modeling  framework  by  ex¬ 
pressing  each  layer’s  functionality  in  probabilistic  terms,  i.e., 
in  the  PHY  layer  we  are  concerned  with  the  probability  of  a 
successful  frame  transmission,  whereas  in  the  MAC  layer  we 
are  concerned  with  its  scheduling  rate  (transmission  proba¬ 
bility).  We  model  the  interaction  between  a  node  and  every 
other  node  in  the  network  by  means  of  interference  matrices. 

3.1  Impact  of  the  PHY  Layer: 

Successful  Lrame  Reception  Probability 

The  quality  of  the  radio  link  between  any  two  nodes  in 
an  ad  hoc  network  may  vary  considerably  over  time,  even 
in  the  case  when  the  two  nodes  are  within  range  of  each 
other.  The  time-variant  nature  of  a  radio  link  is  due  to 
physical  layer  aspects,  as  well  as  the  transmission  activity 
of  all  other  nodes  in  the  network,  which  is  controlled  by 
the  MAC  scheme.  Individually,  a  node’s  transmission  may 
not  interfere  with  some  packet  reception  at  a  distant  node. 
However,  collectively,  the  aggregate  interference  from  many 
simultaneous  transmissions  might  do. 

With  few  exceptions  (e.g.,  [48,  49]),  most  prior  models  of 
MAC  protocols  assume  that  all  packets  received  at  the  same 
time  by  a  receiver  “collide”  and,  therefore,  are  destroyed. 
This  assumption  is  reasonable  if  the  powers  of  the  received 
signals  are  nearly  the  same.  In  most  situations,  however, 
the  powers  of  the  received  signals  are  subject  to  large-scale 
path  loss  propagation,  shadowing,  and  small-scale  multipath 
fading  [50] .  Even  when  packets  from  different  nodes  overlap 
in  time,  it  may  still  be  possible  to  successfully  decode  the 
packet  with  the  strongest  received  signal  strength — the  so- 
called  capture  phenomenon  [51].  Moreover,  a  number  of 
choices  in  system  design  at  the  PHY  layer  have  a  direct 
impact  on  the  successful  reception  of  a  packet  by  a  given 
receiver. 

For  all  the  above  considerations,  a  link  in  a  wireless  ad  hoc 


network  may  generally  be  seen  as  an  asymmetric  channel, 
and  the  conditions  for  a  packet  to  be  successfully  received 
at  one  end  of  the  link  may  be  significantly  different  at  the 
other  end.  Therefore,  we  need  to  consider  network  topology 
in  a  broader  sense,  in  which  asymmetries,  unpredictability, 
and  physical-layer  aspects  are  explicitly  incorporated  into 
the  model. 

For  the  purposes  of  this  paper,  we  classify  the  channel  sig¬ 
naling  methods  into  two  types  [52]:  narrowband  and  spread- 
spectrum  systems.  In  narrowband  signaling,  data  bits  are 
modulated  directly  onto  the  carrier.  On  spread-spectrum 
systems,  some  form  of  coding  is  used  to  get  a  much  wider 
bandwidth  than  narrowband  schemes  do  for  the  same  data 
rate.  In  many  cases,  spread-spectrum  signaling  is  the  signal¬ 
ing  method  of  choice  because  of  its  anti-jamming  capabili¬ 
ties,  robustness  to  multi-path  effects,  lower  power  spectrum 
density,  and  potential  for  multiuser  access  through  CDMA 
techniques,  i.e. ,  the  use  of  orthogonal  codes  that  can  overlap 
in  time  with  little  or  no  effect  on  each  other.  As  an  exam¬ 
ple,  the  IEEE  802.11  standard  supports  both  direct  sequence 
spread  spectrum  (DSSS)  and  frequency-hopped  spread  spec¬ 
trum  (FHSS)  systems  [53]. 

In  this  paper,  we  assume  a  direct  sequence  spread  spec¬ 
trum  ad  hoc  network.  In  such  networks,  multiple  access  in¬ 
terference  (MAI)  plays  a  major  role  on  network  capacity  [37], 
[54]:  in  the  demodulation  of  each  signal,  signals  from  other 
nodes  transmitting  simultaneously  over  the  channel  appear 
as  interference.  The  level  of  such  interference  varies  accord¬ 
ing  to  the  number  of  active  nodes  at  any  time  (dictated  by 
the  underlying  medium  access  control  scheme).  For  this  rea¬ 
son,  not  only  capacity,  but  many  QoS  measures  in  wireless 
ad  hoc  networks  are  intrinsically  dependent  on  the  received 
signal-to-interference-plus-noise  density  ratio  (SINR).  More 
specifically,  let  V  denote  the  finite  set  of  |  |  =  n  nodes 
spanning  the  network  under  consideration,  and  Vr  C  V  the 
subset  of  nodes  whose  signal  powers  can  be  perceived  at 
node  r  significantly.  In  other  words,  by  letting  Vr  C  V  we 
also  consider  the  cases  where  there  exist  nodes  in  V  whose 
contributions  to  the  aggregate  interference  at  node  r  are 
practically  neglible  because  of  the  RF  propagation  effects  in 
place.  Let  now  Pk  denote  the  received  signal  power  at  node 
r  for  a  signal  transmitted  by  node  k  £  Vr.  Then,  the  signal- 
to-interference-plus-noise  density  ratio  SINR(  for  a  signal 
transmitted  by  node  i  and  received  at  node  r,  is  given  by 
(using  a  conventional  matched  filter  receiver)  [55]: 


SINR"  = 


P[  Li 


E  jeVr, +°-r 


(1) 


where  Li  is  the  spreading  gain  (or  bandwidth  expansion  fac¬ 
tor)  of  the  spread-spectrum  system,  is  the  background  or 
thermal  noise  power  at  the  front  end  of  the  receiver  r,  and 
Xj  is  an  on/off  indicator,  i.e., 


{1,  if  j  transmits  at  the  same  time, 
0,  otherwise. 


(2) 


In  Eq.  (1),  it  is  assumed  that  nodes  use  long,  orthogonal 
pseudo-random  sequences  [55],  and  that  the  thermal  and 
background  noise  is  modeled  as  a  white  Gaussian  noise  pro¬ 
cess.  The  same  assumptions  are  used  in  such  simulators  as 
Glomosim,  Ns- 2,  and  Qualnet  [5]. 

Each  time  node  i  transmits  a  frame  to  node  r,  the  (in¬ 
stantaneous)  MAI  level  at  node  r  depends  on  the  nodes  in 
Vr  that  are  transmitting  concurrently  with  node  i,  as  indi¬ 


cated  by  the  variable  X:i ,  which  is  Bernoulli-distributed  with 
probability  Tj,  the  probability  that  node  j  transmits  a  frame 
at  any  time,  according  to  the  MAC  protocol  in  place. 

If  \Vr\  =  nr,  there  are  exactly  2rir~1  combinations  of  ac¬ 
tive  transmitting  nodes  ( interferers )  in  Vr,  excluding  the 
transmitter  i  itself.  In  what  follows,  let  {Cik}  k=1  2«r-i 

denote  the  set  of  such  combinations,  and  C'  be  a  random 
variable  that  indicates  the  occurrence  of  a  specific  combina¬ 
tion  c\k  of  interferers. 

For  simplicity,  we  assume  that,  when  node  i  transmits 
a  frame  to  node  r,  the  set  of  interferers  remain  the  same 
throughout  the  entire  transmission  of  the  frame.  In  reality, 
of  course,  some  nodes  may  become  active  or  inactive  during 
the  course  of  a  frame  transmission.  However,  this  assump¬ 
tion  is  reasonable  if  frames  are  short  and  transmission  rates 
are  high.  Consequently,  during  the  reception  of  a  frame,  and 
for  a  fixed  set  of  interferers,  the  bit-to-bit  variations  that  may 
occur  in  SINR(  will  result  from  RF  propagation  effects  only. 
Typically,  the  received  signal  powers  Pk ,  k  £  14,  are  sub¬ 
ject  to  large-scale  path  loss  ,  shadowing,  and  small-scaling 
multipath  fading  [50]. 

Given  the  considerations  above,  the  probability  q [  that 
a  frame  transmitted  by  i  is  successfully  received  at  r  can 

be  obtained  by  considering  the  set  {c^} k=i . j«r-i  of  all 

possible  combinations  of  active  nodes  in  Vr,  as  follows: 

qi  =  P{  successful  frame  reception  } 

=  J2kp{  successful  frame  reception,^  =  clk } 

=  J2kp{  succ.  frame  reception!  C[  =  crik}P{Cl  =  crik} 

=  J2J(crik)p{C:  =  crik},  (3) 

where  we  wrote  /(clfc)  to  emphasize  the  fact  that  P{succ. 
frame  reception  |  Cf  =  crik}  is  a  function  of  the  specific  com¬ 
bination  c\k  of  interferers.  Its  functional  form  will  depend 
on  the  specific  choice  of  radio  channel  model  and  such  PHY- 
layer  aspects  as  modulation  and  demodulation  schemes,  chan¬ 
nel  coding,  receiver  design,  and  the  like. 

The  above  formalism  allows  for  the  consideration  of  any 
radio  channel  model  and  PHY-layer  aspect  for  computation 
of  the  probability  /(cjj.)  of  successful  frame  reception  con¬ 
ditioned  on  a  certain  MAI  level.  For  instance,  one  could 
choose  to  model  the  impact  of  a  specific  modulation  scheme 
under  correlated  bit  errors  during  frame  reception. 

Now,  let  us  consider  the  probability  that  the  set  of  active 
interferers  is  c)*.,  i.e.,  P{Cf  =  crik\.  This  probability  is  a 
function  of  the  MAC-dependent  transmission  probabilities 
Tj,  j  £  Vr,  which,  for  each  node,  depends  on  the  perceived 
activity  of  other  nodes.  Hence,  assuming  that  Tj,  j  £  Vr, 
is  known,  we  make  the  simplifying  assumption  that  P{C[  = 
elk }  is  given  by 

m=<4}=  II  (i-Tm)  II  T”’  (4) 

n^cik 

where  crik  denotes  the  complement  set  of  crik,  i.e.,  Vr  —  {crik\. 
Finally,  from  Eqs.  (3)  and  (4), 

Qi  =^2,f{cik)  (1  -Tm)  n  Tn ■  (5) 

k  mecVk  n^cik 

3.2  Impact  of  the  MAC  Layer: 

The  Transmission  Probabilities 

So  far,  we  have  assumed  that  a  node  k  £  V  transmits 
a  frame  according  to  some  scheduling  transmission  rate  (or 


transmission  probability)  Tk  dictated  by  the  MAC  protocol. 
In  traditional  analysis  of  MAC  protocols,  the  underlying 
scheduling  point  process  is  usually  assumed  to  be  Poisson- 
distributed,  and  performance  evaluation  of  such  MAC  pro¬ 
tocols  are  usually  conducted  by  allowing  the  offered  load 
to  vary  over  a  certain  range  [6],  [7],  [9],  [8],  [17].  In  these 
models,  it  is  normally  assumed  that  retransmissions  occur 
after  a  randomized  period  in  time  that  is  “sufficiently  long” 
in  order  to  preserve  the  Poisson  assumption  [7].  By  doing 
so,  retransmissions  are  basically  treated  as  new  transmis¬ 
sions,  independent  of  previous  attempts.  Furthermore,  in 
some  cases,  the  packet  length  of  each  retransmission  is  as¬ 
sumed  to  have  been  drawn  from  a  given  distribution,  inde¬ 
pendent  of  the  original  transmission  [52] . 

In  reality,  when  a  reliable  delivery  service  is  desired  at  the 
MAC  level,  a  node  retransmits  a  frame  according  to  some 
specific  rule,  until  the  frame  is  finally  transmitted  or  dis¬ 
carded  after  a  certain  number  of  failures.  At  each  attempt, 
if  we  consider  the  explicit  exchange  of  frames  (ACKs)  ac¬ 
knowledging  the  data  frames  transmitted,  two  events  need 
to  happen  for  a  node  i  to  consider  its  frame  successfully 
transmitted  to  a  node  r:  the  successful  reception  of  i’s  frame 
at  r,  and  the  successful  reception  of  r’s  acknowledgment  at 
i.  Therefore, 

qi  =  P{DATA  successful,  ACK  successful} 

=  P{ACK  succ.  |  DATA  succ.}P{DATA  succ.} 

=  QiQr  =  £*/(<&  )P{C!  =  eik}'Elf&)P{Ctr  =  4l} 

=  Efc  E  J(crik)f(cU)P{Cl  =  clk}P{Ci  =  cU}  (6) 

The  MAC  protocol  dynamically  adjusts  its  behavior  ac¬ 
cording  to  the  feedback  that  is  acquired  during  its  attempts 
to  transmit  a  particular  frame.  The  MAC  protocol  must 
use  the  feedback  information  to  schedule  future  retransmis¬ 
sions  in  a  way  that  minimizes  the  number  of  unsuccessful 
transmissions.  In  this  sense,  MAC  protocols  can  be  seen  as 
dynamic  systems  whose  feedback  information  are  the  suc¬ 
cessful  transmission  probabilities  qt  and  the  corresponding 
outputs  are  the  scheduling  rates  n,  i  £  V.  Therefore,  the 
MAC  operation  can  be  represented  by  some  function  hi(-) 
that  maps  the  feedback  input  into  the  desired  outputs,  which 
are  the  scheduling  rates.  We  use  the  subscript  i  in  the  map¬ 
ping  function  hi(-)  to  denote  a  node-specific  MAC-protocol 
output. 

Strictly  speaking,  the  MAC  operation  should  be  seen  as 
a  time-variant  system,  where  probabilities  should  also  be 
time- variant  (a  stochastic  process).  However,  we  will  assume 
that,  in  the  long-run  (or  in  steady-state),  the  operation  of 
the  MAC  protocol  can  be  represented  by  a  time-invariant 
(linear  or  non-linear)  function  hi(-)  relating  the  probabilities 
q-i  and  r»,  *  €  V,  i.e., 

n  =  hi(qi),  i  €  V.  (7) 

Notice  that,  rigorously  speaking,  rt  may  be  a  function  of 
not  only  qi,  but  also  of  the  steady-state  probability  that  a 
packet  is  available  at  a  node’s  queue,  i.e.,  that  the  queue 
is  not  empty.  Consequently,  n  should  also  be  a  function  of 
the  particular  input  data  traffic  distribution  at  node  i  €  V. 
In  the  literature,  however,  the  throughput  of  a  packet-radio 
network  is  studied  by  assuming  saturation  at  all  nodes,  i.e., 
all  buffers  are  full  [52].  In  this  paper,  we  also  assume  that 
a  packet  is  always  available  for  transmission  at  the  head  of 
a  node’s  queue  (being  a  retransmission  or  not).  For  this 


reason,  we  can  adopt  Eq.  (7)  to  represent  the  steady-state 
operation  of  the  MAC  protocol. 

The  coupled  nonlinear  multivariate  system  on  the  vari¬ 
ables  n  and  qi,  i  £  V ,  represented  by  Eqs.  (4),  (6),  and  (7), 
allows  us  to  consider  the  interdependence  between  the  PHY 
layer  and  the  MAC  layer  explicitly,  as  well  as  the  interde¬ 
pendence  among  all  nodes  in  the  network!  However,  solving 
such  a  nonlinear  system  in  arbitrary  topologies  constitutes 
an  extremely  complex  task  if  we  consider  the  number  of  pos¬ 
sibilities  for  the  functional  forms  of  /(•),  hi(-),  and  network 
sizes.  Furthermore,  checking  for  the  existence  of  a  solution 
for  the  system  can  be  very  cumbersome,  if  not  impossible. 
Hence,  to  simplify  our  modeling  problem  and  to  better  un¬ 
derstand  the  effects  of  different  protocol  parameters  on  the 
probabilities  n  and  qi,  in  the  next  section  we  linearize  the 
above  nonlinear  system. 

3.3  The  Impact  of  Topology: 

The  Linear  Model 

Let  us  first  consider  the  probability  n  that  a  node  trans¬ 
mits  a  packet  at  any  given  time  as  expressed  by  Eq.  (7).  Ac¬ 
cording  to  our  previous  remarks,  the  MAC-dependent  func¬ 
tional  form  hi(-)  may  be  a  linear  or  non-linear  function  of  the 
successful  transmission  probability  qi.  To  generalize,  let  us 
assume  that  it  is  a  non-linear  function.  Therefore,  if  hi(qi)  is 
a  function  with  a  continuous  nth  derivative  throughout  the 
interval  [0,  1],  the  Taylor  series  expansion  of  hi(qi)  is  given 
by 

„n  —  1 

hi(qi)  =  hi( 0)  +  qih'i( 0)  +  •  •  •  +  ^  jyy ^4"  ^(0)  +  Rn, 

(8) 

where  Rn  =  ^  h\n\e),  0  <  e  <  qi. 

Let  hi(qi)  be  a  function  such  that  hi( 0)  «  0,  i.e.,  the  prob¬ 
ability  that  node  i  attempts  to  transmit  a  packet  is  nearly 
zero  when,  in  the  long-run,  qi  =  0.  Such  situation  arises,  for 
instance,  in  the  reasonable  case  where  the  MAC  behavior 
(represented  by  the  function  h(-))  is  such  that  it  does  not 
schedule  transmissions  (i.e.,  its  transmission  probability  is 
almost  zero)  if  it  knows  that  packets  will  not  be  delivered 
(from  the  “feedback  information”  qi  =  0).  In  this  case,  a 
first-order  approximation  for  hi(qi)  is  simply 

Ti  =  hi(qt)  «  aqt,  where  a  = /i'( 0),  (9) 

which,  considering  all  nodes  in  the  topology,  can  be  rewrit¬ 
ten  in  matrix  notation  as  follows: 

r  =  aq,  (10) 

where  r  =  [n  r2  ...  rn]T  and  q  =  [qi  q2  ...  qn]T ■ 

Let  us  now  consider  the  probability  qi  that  a  successful 
transmission  occurs.  According  to  Eq.  (6),  qi  depends  on  a 
number  of  terms  involving  the  transmission  probabilities  r, , 
as  given  by  Eq.  (4).  When  contention  among  nodes  is  high, 
Tj  is  expected  to  be  generally  small.  In  addition,  because 
the  conditional  probability  of  successful  frame  transmission 
/(•)  decreases  as  the  MAI  level  increases  (the  SINR  deteri¬ 
orates),  the  term  in  the  double  summation  of  Eq.  (6)  that 
contributes  the  most  to  the  computation  of  qi  is  exactly  the 
one  that  corresponds  to  the  case  when  none  of  the  potential 
interferers  of  nodes  i  and  r  transmit,  i.e.,  the  MAI  is  null.  In 
this  case,  all  terms  in  the  products  of  Eq.  (4)  will  be  of  the 
form  (1  —  Tm),  and  the  SINR  will  be  the  highest  possible, 


maximizing  /(•).  If  we  let  cl0  and  c(.0  denote  the  combina¬ 
tions  corresponding  to  the  case  when  no  interferer  of  both  r 
and  i  transmit,  i.e.,  cri0  =  clr0  =  {0},  meaning  that  cri0  =  Vr 
and  clr 0  =  Vi,  then  we  can  approximate  qi  as  follows: 

«  =  EE  f(crik)f(cU)P{Cl  =  crik}P{Ci  =  cU} 

k  l 

«  f(cri0)f(ci.0)P{C;  =  cl 0}P{Cl  =  etc} 

=  f(cio)f(cr0)  n  (i  -  Tj)  n  (i  -  Tk),  (ii) 

j£Vr  k£Vi 

From  Eq.  (9), 

Qi  =  /(<4)/(4o)  n(!“  aqj )  Ha"  aqk ^  (12) 

jevr  keVi 

If  we  assume  a  «  1,  and  because  0  <  qi  <  1,  we  can 
approximate  the  previous  products  as  follows: 


Qi  «  f(clo)f(cr0)  |  1  -  a  y  Qi  )  (  1  -  «  E  qk  ] 
V  jevr  )  \  keVi  ) 


f(clo)f{Cro)  1  -  aqA  ■ 

V  je  VrUVi  I 


(13) 


Now,  making  / (c”0) f (c(,0)  =  n»j  and  considering  the  set 
V  of  all  nodes  in  the  network,  with  n  denoting  the  intended 
receiver  of  node  i,  we  have  that 


Ql  =  7T1  -  U7T i  y  Qj  (14) 

lev riuvi 

<?2  =  7T2  -  07T2  Qj  (15) 

ievr2uv2 

Qn  =  nn  —  ann  gy.  (16) 

jevrn  uv„ 

which,  in  matrix  notation  becomes 


Qi  ' 

'  TTi  ‘ 

0 

012 

013  ' 

'  0177 

"  Qi  " 

Q2 

7T2 

021 

0 

023  ' 

0277 

Q2 

Q3 

= 

7T3 

- 

031 

032 

0  ■ 

0377 

Q3 

Qn 

TTn 

071 1 

0772 

0713  • 

0 

Qn 

(17) 

where 

J  am,  if  j  €  Vi  U  Vn 
***  ~  1  0,  if  j  i  Hi  U  Vu 

(18) 

Or,  in  a  more  compact  form, 

q  =  7T  —  <Fq, 

(19) 

leading  to  the  linear  system 

(I  +  f>)q  =  7T, 

(20) 

The  linear  system  in  Eq.  (20)  has  a  solution  if  and  only  if, 
the  vector  it  is  in  the  column  space  of  the  matrix  \F  =  I  +  <F, 
i.e.,  it  is  a  linear  combination  of  the  columns  of  \F.  Given 
the  generality  of  the  matrix  \F ,  which  can  change  every  time 
we  consider  a  different  network  topology  (through  the  inter¬ 
ference  matrix  $),  we  need  to  find  out  which  matrices  IF 
allow  the  linear  system  of  Eq.  (20)  to  have  a  solution.  In 
fact,  we  will  aim  higher  and  ask,  ultimately,  if  the  system  of 
Eq.  (20)  has  a  solution  regardless  of  network  topology.  In¬ 
deed,  the  answer  to  this  question  is  yes  and  is  the  result  of 
Theorem  1  below. 

Definition  1.  A  square  matrix  A  is  diagonally  domi¬ 
nant  if  the  absolute  value  of  each  diagonal  element  is  greater 
than  the  sum  of  absolute  values  of  the  non-diagonal  elements 
in  its  row.  That  is, 

J2  Hbi)l  <  HM)I, 

3 

Property  1.  If  the  diagonal  elements  of  a  square  ma¬ 
trix  A  are  all  positive  and  if  A  and  Ar  are  both  diagonally 
dominant,  then  A  is  positive  definite. 

Theorem  1.  Given  n  >  1,  if  a  <  (n  —  1)_1  the  matrix 
^F  =  I  +  $  is  nonsingular  regardless  of  network  topology. 

Proof.  By  construction,  the  diagonal  elements  of  \F  are 
all  positive  (all  ones).  In  this  case,  for  the  matrices  VF  and 
\Ft  to  be  diagonally  dominant,  we  must  have 

<  1,  for  alH  and  (j>ij  <  1,  for  all  ji,  (21) 

3  i 

where  < pij  are  defined  according  to  Eq.  (18).  Now,  let  K 
denote  the  set  of  non-zero  off-diagonal  elements  in  a  specific 
line  or  column  of  \F.  Because  0  <  70  <  1,  every  line  or 
column  of  \F  satisfies 

y,  ank  <  a7Tmax  <  a  <  (n  —  l)a. 

keK  keK  k£K 

From  the  conditions  in  the  inequalities  in  (21),  and  by  Prop¬ 
erty  1,  it  follows  that  \F  is  positive  definite  if  a  <  ( n —  l)-1. 
Consequently,  by  the  positive  definiteness  property,  \F_1  ex¬ 
ists  and  the  theorem  is  true.  □ 

Considering  the  fact  that  a  is  a  function  of  the  parameters 
of  the  underlying  MAC  protocol  (from  the  linearization  of 
hi(qi)),  the  condition  stated  in  Theorem  1  actually  suggests 
a  design  rule  for  choosing  “feasible”  MAC  parameters  for  a 
given  network  of  size  n.  Note  that,  the  condition  imposed 
on  a  in  Theorem  1  refers  to  the  extreme  case  when  all  nodes 
in  a  network  of  size  n  can  hear  each  other  and  packets  are 
not  missed  due  to  physical  layer  aspects  (ni  =  1  Vi),  i.e., 
a  fully-connected  network  under  ideal  channel  conditions. 
Finally,  if  we  substitute  Eq.  (20)  into  Eq.  (10)  we  get  the 
transmission  probability  vector  r,  given  by 

r  =  a(  I  +  3>)-17r.  (22) 


where  I  is  the  identity  matrix,  $  is  the  n  x  n  matrix  in 
Eq.  (17),  and  it  =  [vri  7r2  . . .  ivn]T .  In  its  essence,  the  matrix 
conveys  all  the  information  regarding  how  each  node  in¬ 
terferes  with  every  other  node  in  the  network  based  on  the 
effect  of  the  PHY  and  MAC  layers.  For  this  reason,  we  refer 
to  fF  as  the  interference  matrix. 


Observe  that,  in  the  worst  case  scenario,  the  linear  sys¬ 
tem  of  Eq.  (20)  can  be  computed  in  up  to  0(n3)  operations 
(where  n  is  the  number  of  nodes).  Depending  on  the  symme¬ 
try  of  the  problem,  the  complexity  can  be  reduced  to  0(n) 
[56].  This  indicates  that  our  modeling  framework  is  quite 
scalable  with  the  number  of  nodes. 


4.  APPLICATION: 

MODELING  THE  IEEE  802.11  DCF 

To  show  how  our  modeling  framework  can  be  applied  to 
the  study  of  MAC  protocols  operating  in  ad  hoc  networks, 
we  apply  our  model  to  the  analysis  of  ad  hoc  networks  using 
the  IEEE  802.11  distributed  coordination  function  (DCF). 

Recently,  we  have  derived  closed-form  expressions  for  the 
first  two  moments  of  a  node’s  service  time  in  a  saturated 
IEEE  802.11  DCF  network,  as  a  function  of  the  channel 
state  as  perceived  by  each  node  [3] .  In  this  model,  the  state 
of  the  channel  is  conveyed  in  the  form  of  channel  state  prob¬ 
abilities,  and  is  expressed  in  terms  of  three  mutually  ex¬ 
clusive  events:  Ei  =  (idle  channel},  Ec  =  {collision},  and 
Es  =  {successful  transmission},  which  are  the  three  events 
that  dominate  the  behavior  of  the  binary  exponential  backoff 
algorithm  in  the  IEEE  802.11  DCF.  To  compute  the  channel 
state  probabilities,  we  use  the  work  by  Bianchi  [4] .  Bianchi’s 
work  provides  a  nonlinear  system  of  equations  relating  the 
steady-state  probabilities  that  a  node  transmits  a  packet  at 
any  given  time,  and  the  respective  collision  probability  of 
this  packet  for  the  case  of  a  fully -connected  network  under 
ideal  channel  conditions  (i.e.,  no  hidden  terminals,  capture, 
or  considerations  of  the  physical  layer  whatsoever). 

In  this  section,  we  extend  Bianchi’s  model  [4]  and  our 
own  prior  model  aimed  at  fully-connected  networks  [3]  to 
consider  a  general  multihop  ad  hoc  network  using  802.11 
DCF  under  non-ideal  channel  conditions.  For  this  purpose, 
we  first  review  the  expressions  for  the  average  service  time 
and  throughput  for  a  saturated  ad  hoc  network  we  derive 
in  [3].  For  the  average  backoff  time,  we  found  that 
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where 
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q  -  2m(l  -  q)m+1 
1-2(1  -q) 


(24) 


In  the  equations  above,  Wmin  is  the  minimum  contention 
window  size  specified  for  the  backoff  operation,  m  is  the 
standard-defined  maximum  power  used  to  set  up  the  maxi¬ 
mum  contention  window  size,  i.e.,  IFmax  =  2mWmin,  q  is  the 
conditional  probability  of  successful  handshake  (assumed 
constant),  and  a  =  api  +tcpc  +  tsps,  where  p;  =  [p]  pzc  p\]T 
is  the  vector  of  channel  state  probabilities  for  the  events 
that  node  i  can  perceive  during  its  backoff  operation,  with 
a,  ts,  and  tc  being  their  corresponding  average  time  dura¬ 
tion.  Given  the  backoff  time  characterization,  the  average 
service  time  equals 


T  =  Tb  +  Ts, 


(25) 


where  Ts  is  the  average  time  to  successfully  transmit  a  packet 
at  the  end  of  the  backoff  operation  (dependent  on  the  packet 
size).  Finally,  since  the  network  is  assumed  to  be  saturated, 
the  average  throughput  per  node  is  simply  1 

s  _  E{Data  Payload}  ^ 

One  of  the  drawbacks  in  our  model  for  fully-connected 
networks  [3]  is  the  fact  that  frames  are  allowed  to  back¬ 
off  infinitely  in  time,  as  opposed  to  what  is  defined  in  the 
standard,  where  retransmission  counters  limit  the  number 

1For  more  details,  see  [3]. 


of  attempts  to  transmit  a  particular  frame,  after  which  the 
frame  is  dropped.  The  infinite  backoff  abstraction  makes 
the  model  much  more  tractable,  but  makes  it  too  conser¬ 
vative,  predicting  higher  service  times  (and  consequently, 
lower  throughputs)  than  what  actually  happens  in  the  IEEE 
802.11  DCF  [3].  Accordingly,  we  extend  our  prior  work  and 
develop  the  average  service  time  for  the  case  of  a  finite  back¬ 
off  operation. 

We  keep  our  assumption  of  a  constant  conditional  success¬ 
ful  handshake  probability  q  per  transmission  attempt,  i.e., 
at  the  end  of  a  backoff  stage,  the  probability  of  a  successful 
handshake  is  constant,  regardless  of  the  number  of  previous 
attempts.  We  should  point  out  that,  later  in  our  develop¬ 
ment,  this  constant  probability  will  depend  on  the  specific 
receiver,  as  expressed  in  Eq.  (6). 

Now,  let  M  be  the  maximum  number  of  times  a  frame 
can  be  retransmitted,  i.e.,  the  maximum  number  of  backoff 
stages  a  frame  can  undergo.  In  this  case,  we  have  a  truncated 
geometric  distribution,  given  by 

P{B  =  k}=  k  =  1,2,  ... .  ,  M,  (27) 

Hence,  following  our  development  in  [3],  the  average  ser¬ 
vice  time  is  now  given  by  (omitting  intermediate  steps,  which 
are  straigthforward  but  tedious): 

Tb  =  +  03  tc,  (28) 

Ai  +  A2  +  A3 
1  -  (1  -  q)M 

2g{i  —  [2(1  —  g)]m}  _1  +  (1_a)m 

2q  -  1  1  q> 

(2m+1  -  1)  (1  -  q )m  [l  -  (1  -  q)M-mj 

2m  { (1  —  9)m+1  -  (1  -  g)M[l  +  q(M  -m-  1)]} 
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i  -  (i  -  g)M(i  +  qM) 

g[l-(l-g)M] 

(l_g)_(l_g)M[l  +  g(M_l)] 

g[i —  (i  —  g)M] 

where  the  parameters  a,  m,  Wmin,  and  tc  are  defined  as 
before. 

To  find  the  channel  probabilities  for  a  general  multihop  ad 
hoc  network,  we  use  the  work  by  Bianchi  [4]  as  our  starting 
point.  Bianchi  provided  a  model  to  evaluate  the  saturation 
throughput  of  the  IEEE  802.11  DCF  MAC  protocol  under 
the  hypothesis  of  ideal  channel  conditions,  i.e.,  no  hidden 
terminals  and  capture.  The  novelty  of  his  approach  comes 
from  the  fact  that  many  of  the  details  of  the  IEEE  802.11  bi¬ 
nary  exponential  backoff  algorithm  were  taken  into  account 
when  modeling  the  stochastic  process  describing  the  back¬ 
off  time  counter  of  a  given  node.  To  accomplish  this,  the 
key  assumption  adopted  in  the  model  was  that  each  frame 
collided  with  a  constant  and  independent  probability  p  at 
each  transmission  attempt,  regardless  of  the  number  of  re¬ 
transmissions  already  occurred.  This  probability  was  called 
the  conditional  collision  probability,  i.e.,  the  probability  of 
a  collision  experienced  by  a  frame  being  transmitted  in  the 
channel.  In  Bianchi’s  development,  a  fixed  number  n  of 
nodes  was  assumed,  with  each  node  always  having  a  packet 
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available  for  transmission,  i.e.,  the  transmission  queue  of 
each  node  was  always  nonempty.  Then,  the  probability  r 
that  a  node  transmits  a  frame  at  any  time  was  obtained 
by  modeling  the  stochastic  process  representing  the  back¬ 
off  time  counter  as  a  bidimensional  discrete-time  Markov 
process.  For  convenience,  we  repeat  the  expression  for  r 
obtained  by  Bianchi: 


_ 2(1  ~  2 p) _ 

(1  -  2p)(Wmin  +  1)  +  pWmin(l  -  (2 p)m)  ' 


(29) 


at  pi  =  0  is  given  by 
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+  0(Pi), 


(32) 


where  0(p2)  accounts  for  the  second  and  high  order  terms 
in  the  Taylor  series  expansion.  Hence,  a  first  order  approx¬ 
imation  of  r/i(pi)  is  simply 


ViiPi) 


Wmin  -  1  2  Wmin 

Wmin  +  1  (Wmin  +  l)2 


(33) 


To  model  a  multihop  wireless  ad  hoc  network  using  the 
IEEE  802.11  DCF,  we  need  to  consider  each  node  individ¬ 
ually.  As  before,  let  V  denote  the  finite  set  of  \V\  =  n 
nodes  spanning  the  network  under  consideration.  Each  time 
a  node  i  £  V  attempts  to  perform  a  handshake  (at  the  end  of 
a  backoff  stage)  it  will  experience  a  failure  with  probability 
Pi.  As  in  Bianchi’s  work,  and  given  that  we  are  not  consid¬ 
ering  mobility,  we  will  also  assume  that  pi  is  constant  and 
independent  of  the  number  of  retransmissions  that  already 
occurred  (constant  with  respect  to  the  number  of  retrans¬ 
missions,  but  dependent  on  the  chosen  receiver  because  we 
are  including  the  physical  layer  aspects  now).  In  the  specific 
case  of  single-hop  networks  under  ideal  channel  conditions, 
previous  models  implicitly  assumed  that  a  successful  recep¬ 
tion  of  an  RTS  (DATA)  frame  was  sufficient  for  a  successful 
handshake  to  occur.  However,  in  the  multihop  environment 
under  non-ideal  conditions,  a  successful  handshake  depends 
on  the  occurrence  of  at  least  two  events:  the  successful  recep¬ 
tion  of  both  RTS  (DATA)  and  CTS  (ACK)  frames.  As  far 
as  the  backoff  dynamics  is  concerned,  only  the  reception  of 
RTS  and  CTS  packets  determines  if  a  node  leaves  the  back¬ 
off  state  or  not.  For  the  moment,  we  will  disregard  possible 
collisions  involving  DATA  packets  or  ACK  control  packets 
in  the  four-way  handshake  scheme.  For  this  reason,  we  will 
refer  to  the  probability  pi  in  a  more  general  sense,  as  the 
conditional  probability  of  a  failed  handshake.  Given  pi  and 
the  previous  assumptions,  the  stochastic  process  governing 
a  node’s  backoff  time  counter  is  well  defined  and,  therefore, 
its  steady  state  transmission  probability  n  can  be  expressed 
in  the  same  way  as  in  Eq.  (29),  except  that  it  now  refers  to 
each  node  individually: 

=  _ 2(1  -  2pQ _ 

1  (1  —  2pi)(Wmin  +  1)  +PiWmin(l  —  (2pi)m)'  1  ’ 

Eq.  (30)  is  exactly  what  we  need  to  apply  our  model¬ 
ing  methodology:  it  gives  us  the  functional  form  hi(qi)  by 
which  the  MAC  layer  relates  the  steady-state  transmission 
probability  n  with  the  successful  transmission  probability 
q.i  =  1  —  Pi.  Therefore,  all  that  is  left  for  us  to  do  is  to  com¬ 
pute  a  first  order  approximation  for  it.  For  this  purpose, 
let  Tji  be  the  probability  that  node  i  does  not  transmit  in  a 
randomly  chosen  slot  time,  i.e., 

(1  -  2Pi)(Wmin  -  1)  +  PilFminU  -  (2pi)m) 

Vi  ^  (1  -  2 PiXWmin  +  1)  +  Pi  ITmin  (1  -  (2pi)™)  ' 

(31) 

Given  the  continuity  of  both  pi(pi)  and  its  derivatives2  in 
the  interval  pi  £  (0, 1),  the  Taylor  series  expansion  of  ??i(pi) 

2Continuity  with  respect  to  the  critical  value  p  =  1/2  can 
be  shown  by  simply  rewriting  »p(pi)  in  the  same  way  as  it 
was  done  for  n(pi)  by  Bianchi  [4]. 


which,  in  terms  of  qi  =  1  —  Pi  becomes 


Vi{Qi) 


2Wmin  IFmin  +  2Wmin  —  1 

(Wmin  +  l)2  ^  (Wmin  +  l)2 
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Given  that  n  =  1  —  rp,  we  can  finally  express  n  in  terms  of 
the  probability  qt  as  follows: 


n(qi) 


2Wmin 

(Wmin  +  l)2^' 


(35) 


If  we  consider  all  nodes  in  the  topology,  Eq.  (35)  can  be 
rewritten  in  matrix  notation  as  follows: 


r  =  aq,  (36) 

where  r  =  [n  r2  ...  rn]T,  a  =  2Wmm/(Wmm  +  l)2,  and 
q  =  [<?i  12  •••  qn]T. 

We  can  now  deal  with  the  computation  of  the  channel- 
state  probability  vector  pi  =  [p]  pzc  p\]T  for  every  node  i  £ 
V.  For  this  purpose,  we  need  to  consider  the  carrier  sensing 
range  during  the  times  a  node  is  sensing  the  channel.  In 
general,  the  carrier  sensing  range  of  a  node  is  larger  than 
its  transmission  range.  Therefore,  because  the  dynamics 
of  a  node  in  backoff  depends  on  its  perceived  state  of  the 
channel,  we  need  to  take  into  account  the  carrier  sensing 
range  as  opposed  to  the  transmission  range  when  computing 
the  channel  state  probability  vector  p;.  That  said,  let  Ri 
denote  the  set  of  nodes  within  the  carrier  sensing  range  of 
node  i.  At  any  given  time,  the  probability  that  there  exists 
some  node  from  Ri  transmitting  a  frame  while  node  i  is  in 
backoff  is 

Plr  -  1  -  n  0-  -  Ti)  (37) 

jeRi 


The  probability  p\uc  that  a  transmission  within  i’ s  sensing 
range  is  successful  is  given  the  probability  that  some  node 
in  Ri  transmits  successfully,  conditioned  on  the  fact  that  at 
least  one  node  in  Ri  attempted  to  transmits,  i.e., 
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Ptr 
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Xifcefli  IfcTfc 
Ptr 


From  Eqs.  (37)  and  (38),  the  probability  p\  that  a  suc¬ 
cessful  transmission  occurs  within  i’s  sensing  range  is  then 
pi  =  P{ES}  =  p\rplsuc.  Accordingly,  the  probability  p\ 
that  the  channel  within  i’s  sensing  range  is  idle  is  simply 
p\  =  P{Ei}  =  1  —  Ptr 5  an<i  the  probability  pzc  that  a  collision 


occurs  within  i' s  sensing  range  is  given  by  pzc  =  P{EC}  = 

Ptr(.I  Psuc )■ 

Finally,  we  need  to  know  the  values  of  the  time  intervals 
ts  and  tc  which,  for  the  four-way  handshake  mechanism,  are 
given  by 

ts  =  RTS  +  SIFS  +  5  +  CTS  +  SIFS  +  S  +  H  +  E{P} 

+  SIFS  +  5  +  ACK  +  DIFS  +  <5, 

tc  =  RTS  +  DIFS  +  S.  (39) 

where  E{P}  =  P  for  fixed  packet  sizes.  Notice  that  the 
value  of  Ts  in  Eq.  (25),  is  simply  ts  —  DIFS. 

5.  MODEL  VALIDATION 

In  this  Section,  we  evaluate  the  accuracy  of  our  model 
in  predicting  the  performance  of  nodes  in  multihop  ad  hoc 
networks  using  the  IEEE  802.11  DCF,  focusing  on  ad  hoc 
networks  with  static  topologies  under  saturation.  Analytical 
results  are  compared  to  numerical  simulations  using  Qualnet 
v3.5  [5]. 

5.1  Scenarios  Used  for  Comparison 

5.1.1  Radio  Channel  Model 

For  the  purposes  of  validation  of  our  modeling  framework, 
we  consider  a  radio  channel  model  with  large-scale  path  loss 
propagation  effects  only.  We  choose  a  simpler  model  in  order 
to  focus  on  the  model’s  ability  to  faithfully  represent  interde¬ 
pendencies  among  the  nodes  and  the  per-node  performance 
metrics.  Nonetheless,  as  mentioned  in  Section  3.1,  other 
effects  such  as  shadowing  and  small-scale  multipath  fading 
may  be  equally  taken  into  consideration  when  computing 
the  conditional  probabilities  /(c(j.)  of  a  successful  frame  re¬ 
ception  in  Eq.  (3).  The  path  loss  propagation  model  we 
choose  is  the  two-ray  ground  reflection  model  [50]. 

5.1.2  IEEE  802. 1 1  Parameters 

Regarding  the  IEEE  802.11  PHY  layer,  we  use  direct 
sequence  spread  spectrum  (DSSS)  with  a  raw  bit  rate  of 
1  Mbps  with  DBPSK  modulation.  Under  this  configura¬ 
tion,  according  to  the  standard,  both  preamble  and  MAC 
protocol  data  unit  (MPDU)  are  transmitted  at  the  same 
basic  rate,  under  the  same  modulation  scheme  [53].  With¬ 
out  loss  of  generality,  this  is  done  just  to  ease  the  analysis 
and  not  adding  more  complexity  to  the  analytical  model  by 
treating  each  segment  separately,  with  different  modulation 
schemes. 

Each  node  has  the  same  transmit  power  and,  for  the  given 
path  loss  propagation  model,  we  select  receive  and  carrier 
sensing  thresholds  in  a  way  that  the  radio  range  is  set  to 
200  m  and  carrier  sensing  range  is  set  to  400  m.  Table  5.1.2 
summarizes  the  rest  of  the  parameters  used  for  PHY  and 
MAC  layers. 

Because  multipath  fading  is  not  considered,  and  MAI  is 
disregarded  during  the  approximation  of  the  successful  frame 
reception  probability  in  Eq.  (11),  bit-to-bit  error  dependen¬ 
cies  are  practically  null.  Hence,  the  only  kind  of  interfer¬ 
ence  that  is  left  is  the  thermal  and  background  noise,  which 
is  assumed  to  be  white  Gaussian  noise.  Consequently,  we 
can  treat  bit  errors  independently,  as  it  is  assumed  in  Qual¬ 
net  as  well.  Let  7(0(0)  denote  the  SINR  at  node  r  for  a  bit 
transmitted  by  i  when  none  of  r’s  interferers  transmits,  i.e., 
the  combination  c[0,  as  defined  in  Section  3.3.  If  K  is  the 


Table  1:  IEEE  802.11  Simulation  Parameters. 


MAC 

PHY 

lTmin 

32 

Temperature  (Kelvin) 

290 

^max 

1024 

Noise  factor 

10 

MAC  Header  (bytes) 

34 

Transmission  power  (dBm) 

10 

ACK  (bytes) 

38 

Sensitiviy  of  PHY  (dBm) 

-87.039 

CTS  (bytes) 

38 

Minimum  power  for 

RTS  (bytes) 

44 

received  packet  (dBm) 

-76.067 

Slot  Time  (fisec) 

20 

Packet  reception  model 

BER 

SIFS  (/rsec) 

10 

DIFS  (/rsec) 

50 

length  of  the  physical  layer  convergence  protocol  data  unit 
(PPDU),  and  Pb( 7)  is  the  bit-error  probability  of  DBPSK 
for  a  certain  SINR  level  7,  then  the  conditional  probability 
/(cio)  of  successful  frame  reception  in  Eq.  (11)  is  given  by 

/(cro)  =  {l-A[7(Go)]}K.  (40) 

Accordingly,  the  same  procedure  is  used  for  computation  of 
f(4 0)  m  Eq.  (11). 

5.1.3  Interference  Matrix  Computation 

Given  the  computations  of  /(c(0)  and  /(c(.0)  as  above, 
one  can  finally  build  the  linear  system  of  Eq.  (20)  to  com¬ 
pute  the  transmission  probability  vector  r  in  Eq.  (22).  For 
that,  we  need  to  know  who  is  the  set  of  interferers  for  each 
transmitter/receiver  pair,  as  used  in  Eq.  (13). 

Because  aggregate  MAI  is  disregarded  in  the  linear  model, 
we  select  as  a  potential  interferer  of  a  node  r  any  node  j  £  V 
whose  signal  power,  when  received  at  r,  is  above  the  selected 
carrier  sensing  threshold.  Any  node  not  conforming  to  that, 
we  discard  as  an  interferer  of  r  (we  would  certainly  include 
such  signals  if  MAI  were  considered  and  its  contribution 
to  aggregate  interference  were  significant).  Finally,  given 
the  transmission  probability  vector  r,  one  can  compute  the 
channel  state  probability  vectors  p;  as  in  Section  4  and, 
consequently,  per-node  throughputs  according  to  Eq.  (26). 

5.1.4  Simulation  Setup 

Because  of  the  selected  radio  range,  nodes  are  randomly 
placed  in  an  area  of  1000  x  1000  m,  a  size  big  enough  to 
deploy  a  multihop  network.  The  only  constraint  imposed 
is  that  the  resulting  graph  is  connected.  This  constraint 
is  added  simply  to  make  sure  that  all  nodes  have  at  least 
one  neighbor,  and  that  considerable  channel  contention  and 
hidden-terminal  effects  are  present  in  the  scenarios. 

In  simulations,  each  node  chooses  the  same  neighbor  node 
(i.e.,  a  node  within  its  transmission  range)  for  all  its  trans¬ 
missions  using  the  same  CBR  source  rate.  We  pick  source 
rates  high  enough  to  saturate  all  nodes  in  the  network. 
Packet  sizes  are  fixed  to  1500  bytes  (IP  packet)  and  each 
simulation  run  corresponds  to  5  minutes  of  data  traffic.  We 
repeat  the  experiment  for  50  seeds,  with  each  trial  corre¬ 
sponding  to  a  different  initial  transmission  time  at  each 
node.  Initial  transmission  times  are  randomly  chosen  within 
the  interval  [0,  0.01]  s.  This  is  done  to  allow  the  exponen¬ 
tial  backoff  algorithm  to  be  triggered  at  different  instants 
in  time,  at  each  node,  so  that  different  state  evolutions  are 
taken  into  account  for  the  same  topology. 

5.2  Accuracy  of  Models 

Two  aspects  of  interest  in  our  modeling  approach  are  seal- 
ability  and  per-node  performance.  To  illustrate  such  aspects, 
let  us  consider  two  network  topologies:  one  with  50  nodes 


and  another  with  100  nodes.  Figures  1  and  2  show  their 
respective  topologies.  In  the  figures,  the  numbers  inside  the 
areas  indicate  the  node  ID.  To  provide  an  idea  of  the  den¬ 
sity  of  nodes  in  each  topology,  a  line  between  two  nodes 
indicates  that  they  are  within  200  meters  of  one  another. 
As  it  can  be  observed,  the  created  topologies  provide  both 
highly-dense  as  well  as  poorly-connected  areas.  Needless  to 
say,  interference  at  a  given  receiver  in  the  analytical  model 
and  simulations  is  not  confined  to  only  those  sources  shown 
directly  “connected”  in  the  figures. 


Figure  1:  Network  topology  for  the  50-node  net¬ 
work. 


Figure  2:  Network  topology  for  the  100-node  net¬ 
work. 

Figure  3  shows  the  simulation  and  analytical  results  for 
the  50-node  network.  Figures  4(a)  and  4(b)  show  the  results 
for  the  100-node  network.  In  the  graphs,  the  x  axis  contains 
the  node  ID,  and  the  y  axis  shows  the  respective  through¬ 
put.  As  it  can  be  seen,  the  predicted  performance  correlates 
very  well  with  practically  every  node  in  the  topology  (in 
some  cases,  providing  very  close  results).  In  the  cases  where 
the  model  is  off  by  some  factor,  the  predicted  performance 
follows  the  observed  pattern  in  the  majority  of  the  cases. 

A  more  statistically  significant  result  is  obtained  if  we 


Figure  3:  Throughput:  simulations  versus  analytical 
model  for  the  50-node  network. 


evaluate  the  performance  of  our  modeling  approach  over  a 
number  of  different  topologies.  For  this  purpose,  we  con¬ 
sidered  10  network  topologies  with  100  nodes  each,  all  ran¬ 
domly  generated  as  before.  The  topologies  used  for  compu- 
tating  the  histogram  were  similar  to  that  shown  in  Fig.  2. 
We  compute  the  percentage  prediction  error  with  respect 
to  the  maximum  range  of  throughput  values  observed  in 
simulations  (for  that  particular  topology)  for  each  node  in 
each  randomly  generated  topology.  This  way,  we  weigh 
the  prediction  error  with  respect  to  the  dynamic  range  of 
throughput  values  obtained  in  simulations.  We  obtained 
a  histogram  for  each  topology  by  counting  the  number  of 
nodes  within  a  certain  percentage  prediction  error,  and  we 
then  averaged  the  histograms  over  all  topologies.  Figure  5 
shows  our  results.  As  we  can  see,  the  percentage  predic¬ 
tion  error  is  within  20%  in  about  90%  of  the  nodes,  showing 
how  close  our  analytical  model  is  in  predicting  the  results 
obtained  in  simulations. 

Finally,  we  show  the  application  of  our  model  in  predicting 
the  average  throughput  for  a  given  network  topology.  For 
this  purpose,  we  consider  four  topologies  with  25,  50,  75, 
and  100  nodes.  Figure  6  shows  the  comparison  between  the 
average  throughput  obtained  in  both  simulations  and  ana¬ 
lytical  model.  As  we  can  see,  the  analytical  model  provides 
very  good  results  compared  to  simulation  results.  This  last 
figure  also  shows  how  scalable  is  our  model  in  computing 
throughput  for  different  network  sizes  and  parameters. 

5.3  Modeling  Time 

The  importance  of  the  above  results  and  the  strength  of 
our  analytical  model  become  apparent  when  we  analyze  the 
time  required  to  obtain  the  above  results  through  simula¬ 
tion  and  with  our  analytical  model.  In  Qualnet,  each  run 
of  the  simulation  for  the  100-node  scenario  (consisting  of 
a  5-minute  data  traffic)  takes  about  1,118.4  seconds  (19.69 
minutes)  in  a  Sun  Blade  100  machine  running  Solaris  5.8. 
For  the  50  seeds  needed,  this  corresponds  to  16.41  hours 
of  simulation.  In  this  same  machine,  our  analytical  model, 
implemented  in  Matlab  6.0,  takes  about  0.44  seconds.  This 
corresponds  to  a  time  saving  of  more  than  134,000  times. 
Clearly,  this  is  a  strong  argument  supporting  the  case  for 
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Figure  4:  Throughput:  simulations  versus  analytical 
(b)  Last  50  nodes. 


Figure  5:  Percentage  error  prediction  histogram 
over  10  random  topologies. 


the  need  of  powerful  analytical  models  for  multihop  ad  hoc 
networks,  such  as  the  one  we  are  providing  here. 

5.4  Model  Limitations 

The  simplifications  in  our  application  of  the  modeling 
framework  to  the  analysis  of  ad  hoc  networks  using  the 
IEEE  802.11  DCF  could  explain  some  of  the  disagreements 
between  the  analytical  model  and  simulations. 

In  the  model  we  first  presented  in  [3],  it  is  assumed  that 
periods  of  collisions  experienced  by  colliding  nodes  have  the 
same  duration  as  periods  in  which  the  channel  is  sensed 
busy  by  non-colliding  nodes.  This  assumption  has  a  direct 
impact  on  the  number  of  nodes  contending  for  the  channel 
at  any  given  time.  When  two  or  more  nodes  transmit  at 
about  the  same  time,  if  their  packets  collide,  they  will  wait 
CTS_timeout  seconds  until  they  figure  out  that  a  collision 
has  occurred.  On  the  other  hand,  nodes  that  did  not  trans¬ 
mit  at  that  particular  time  will  sense  a  “clear  channel”  much 
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Figure  6:  Average  throughput  for  different  network 


earlier  before  colliding  nodes  do  and,  therefore,  are  more 
likely  to  be  ready  to  transmit  before  the  colliding  nodes. 
For  this  reason,  our  model  implies  a  more  “aggressive”  net¬ 
work,  where  the  probability  of  having  a  node  transmitting 
at  any  given  time  is  higher  than  in  real  scenarios. 

An  important  limitation  of  the  model  has  to  do  with  the 
impact  of  the  capture  effects.  In  the  linearization  procedure 
of  Section  3,  we  only  considered  the  terms  corresponding  to 
maximum  SINR,  which  were  exactly  those  cases  where  none 
of  the  interfering  nodes  were  transmitting  at  the  same  time 
as  the  node  under  consideration.  In  contrast,  the  capture 
effects  are  present  in  Qualnet  simulations,  which  are  known 
to  affect  throughput  performance  considerably  [51].  Never¬ 
theless,  the  analytical  model  seems  to  capture  the  general 
trend  and  provide  a  striking  correlation,  even  in  a  scenario 
with  100  nodes,  where  the  contention  degree  and  hidden 
terminals  are  significant. 

Another  limitation  of  the  modeling  framework  introduced 


ill  Section  3  is  that  the  interference  matrix  of  a  network  cor¬ 
responds  to  a  snapshot  of  the  transmission  decisions  made 
by  the  network  nodes,  i.e.,  it  assumes  a  constant  choice  of  a 
receiver  by  any  one  sender.  In  practice  a  node  will  choose  to 
transmit  to  different  neighbors  or  groups  of  neighbors.  The 
interference  matrix  needs  to  be  extended  to  take  this  effect 
into  account.3 

6.  ADDITIONAL  APPLICATION: 
EVALUATING  THE  EFFECT  OF 
CARRIER-SENSING  RANGE 

Here  we  illustrate  how  our  modeling  framework  can  be 
applied  to  study  different  aspects  of  PHY  and  MAC  layers 
in  ad  hoc  networks  using  the  IEEE  802.11  DCF.  We  focus 
on  the  impact  of  the  carrier  sensing  range  on  network  per¬ 
formance.  In  the  literature,  the  effect  of  these  parameters 
on  network  performance  has  not  been  treated  appropriately 
and,  more  importantly,  no  analytical  model  to  date  has  been 
able  to  deal  with  this  problem  as  faithfully  as  our  modeling 
framework  can. 

In  reality,  none  of  the  PHY-layer  parameters  can  be  “set” 
or  “adjusted”  beforehand.  Despite  the  fact  that  current 
wireless  cards  permit  setting  up  specific  power  thresholds 
for  both  receive  and  sensing  operations,  the  actual  trans¬ 
mission  and  sensing  ranges  depend  on  many  other  factors, 
such  as  path-loss  propagation  (which,  in  turn,  depends  on 
the  surrounding  environment),  multipath  fading,  transmit 
power,  and  antenna  gains.  However,  given  a  radio  propaga¬ 
tion  model  and  a  specific  PHY  layer,  we  can  study  the  im¬ 
pact  of  the  carrier-sensing  range  using  our  analytical  model. 

In  the  analysis  that  follows,  we  apply  the  same  PHY-  and 
MAC-layer  parameters  used  in  Section  5  under  the  same 
channel  conditions  (no  multipath  fading  or  shadowing).  The 
only  exception  is  the  “Sensitivity  of  PHY,”  which  is  a  “tun¬ 
able”  parameter  to  set  the  sensing  range  in  Qualnet  simula¬ 
tions,  for  a  given  radio  propagation  model  and  PHY  layer. 
In  other  words,  the  transmission  power  of  each  node  is  kept 
fixed  (10  dBm,  in  our  case)  and  the  sensing  range  is  set  by 
changing  the  above  threshold  (this  is  how  one  would  proceed 
in  numerical  simulations,  too).  We  use  the  same  topology 
depicted  in  Figure  2. 

To  study  the  impact  of  the  carrier- sensing  range  on  net¬ 
work  performance,  we  hold  constant  the  transmission  range 
and  vary  the  carrier  sensing  range  of  each  node  to  obtain 
the  corresponding  throughput  performance.  We  investigate 
three  carrier  sensing  ranges:  200  m,  300  m,  and  400  m.  Fig¬ 
ure  7  shows  the  results  for  the  average  aggregate  through¬ 
put  for  both  simulations  and  analytical  model.  As  we  can 
see,  when  the  carrier  sensing  range  increases,  the  average 
throughput  decreases  for  both  simulations  and  analytical 
models.  One  reason  behind  this  phenomenon  is  that,  when 
the  carrier  sensing  range  increases,  the  nodes  spend  more 
time  in  backoff,  freezing  their  backoff  counters  more  fre¬ 
quently  as  a  result  of  being  more  sensitive  to  any  activity 
in  the  channel,  due  to  their  higher  sensing  range.  There¬ 
fore,  the  average  service  time  increases  and,  consequently, 
throughput  decreases. 

A  more  intriguing  behavior  renders  the  results  shown  in 
Figure  7.  The  carrier  sensing  range  is  usually  set  to  twice 
the  transmission  range  to  avoid  some  of  the  hidden-terminal 

3This  effort  is  the  subject  of  a  forthcoming  publication. 


problems,  because  nodes  are  able  to  sense  the  activity  of 
other  nodes  not  in  their  transmission  range  during  their 
backoff.  Therefore,  when  the  carrier-sensing  range  is  set 
to  a  value  equal  to  the  transmission  range  (in  our  case, 
200  m),  the  effects  of  the  hidden  terminal  problem  become 
more  apparent,  because  more  RTS/CTS  collisions  happen. 
In  simulations,  this  results  in  a  high  degree  of  unfairness, 
typical  of  the  IEEE  802.11,  rather  than  causing  an  average 
lower  throughput  among  nodes.  In  fact,  in  simulations,  a 
few  nodes  succeed  in  acquiring  the  channel  when  the  carrier 
sensing  range  is  small  and,  therefore,  attain  high  throughput 
values,  but  other  nodes  have  zero  or  practically  null  through¬ 
put.  Therefore,  the  “average”  result  shown  in  Figure  7  for 
the  200-m  sensing  range  actually  refers  to  high  throughput 
values  (of  a  few  nodes),  averaged  over  all  nodes.  Figure  8 
shows  the  simulation  results  on  a  per-node  basis. 

The  analytical  model  we  designed  does  not  incorporate 
the  fairness  problems  of  the  IEEE  802.11.  Consequently,  due 
to  the  hidden-terminal  problems,  nodes  have  lower  through¬ 
put,  but  still  transmit  something.  Figure  9  shows  the  ana¬ 
lytical  results. 
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Figure  7:  The  impact  of  carrier  sensing  range  on 
average  network  throughput. 


7.  CONCLUSIONS 

We  have  introduced  a  novel  modeling  framework  for  mul¬ 
tihop  ad  hoc  networks  that  explicitly  takes  into  account 
the  impact  of  the  PHY  layer  on  the  operation  and  perfor¬ 
mance  of  the  MAC  layer.  A  key  feature  of  our  model  is  that 
nodes  can  be  modeled  individually,  i.e.,  it  allows  a  per-node 
setup  of  many  layer-specific  parameters.  Moreover,  no  spa¬ 
tial  probability  distribution  or  pre-arrangement  of  nodes  is 
assumed.  Instead,  it  allows  the  modeling  of  specific  network 
topologies  and  radio  channel  models. 

Our  analytical  model  constitutes  a  tool  for  a  fast,  accu¬ 
rate,  and  efficient  evaluation  of  a  node’s  performance  for  any 
system  parameter  value  and  topology  of  an  ad  hoc  network. 
We  have  illustrated  this  using  the  IEEE  802.11  DCF  as  an 
example.  Attempting  to  characterize  a  node’s  performance 
only  by  simulations  would  require  considerable  computa¬ 
tional  effort,  specially  when  scalability  is  one  of  the  main 
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Figure  8:  Average  throughput  for  simulations  when 
the  carrier  sensing  range  is  200m. 


Node  ID 

Figure  9:  Average  throughput  for  the  analytical 
model  when  the  carrier  sensing  range  is  200m. 


concerns.  To  be  able  to  draw  general  conclusions,  simu¬ 
lations  would  have  to  be  run  for  a  very  large  number  of 
network  sizes  and  topologies.  Furthermore,  if  we  are  inter¬ 
ested  in  evaluating  the  impact  of  a  single  parameter  to  a 
node’s  service  time,  for  instance,  simulations  would  have  to 
be  repeated  for  each  parameter  value  of  interest. 

Our  model  was  validated  through  discrete  event  simula¬ 
tions  using  the  popular  Qualnet  tool,  and  the  analytical  re¬ 
sults  show  a  striking  correlation  with  the  results  obtained 
via  simulation.  The  importance  of  the  analytical  model  is 
that  the  time  needed  to  obtain  the  desired  results  takes  a 
very  small  fraction  of  the  time  required  to  obtain  the  same 
results  via  discrete-event  simulations. 

Lastly,  we  presented  the  first  analytical  treatment  of  the 
impact  of  the  carrier  sensing  range  on  network  throughput 
in  ad  hoc  networks  using  the  IEEE  802.11  DCF.  As  we  ob¬ 
served,  average  throughput  decreases  as  the  carrier  sensing 
increases  for  a  fixed  transmission  range.  Moreover,  we  have 
noticed  that  the  unfairness  problems  of  the  IEEE  802.11  be¬ 
come  apparent  when  the  carrier  sensing  range  is  equal  to  the 


transmission  range,  where  a  few  nodes  acquire  the  channel  in 
detriment  of  other  nodes,  whose  throughput  is  deteriorated 
to  almost  null  values  due  to  the  increase  of  collisions  from 
hidden  terminals.  The  analytical  model,  on  the  other  hand, 
produces  results  assuming  a  “fair”  operation  of  the  IEEE 
802.11  DCF  and,  consequently,  higher  throughputs  for  that 
case. 
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